A novel detection method for wheat aging based on the delayed luminescence

Wheat aging plays an important role in assessing storage wheat quality and its subsequent processing purposes. The conventional detection methods for wheat aging are mainly involved in chemical techniques, which are time-consuming as well as waste part of wheat samples for each detection. Although some physical detection methods have obtained gratifying results, it is extremely hard to expand their application fields but to stay in the theory stage. For this reason, a novel nondestructive detection model for wheat aging based on the delayed luminescence (DL) has been proposed in this paper. Specifically, after collecting enough sample data, we first took advantage of certain hyperbolic function to fit DL signal, and then used four parameters of the hyperbolic function to feature the decay trend of the DL signal. Secondly, in order to better feature the DL signal, we extracted other six features together with above four features to form the input feature vector. Finally, as the bidirectional long short-term memory (Bi-LSTM) network lacked error-correcting performance, the Bi-LSTM network based on Walsh coding (Walsh-Bi-LSTM) mechanism was proposed to establish the detection model, which made the detection model have the error-correcting performance by reasonably splitting the multi-classification target task. Shown by experimental results, the newly proposed wheat aging detection model is able to achieve 94.00% accuracy in the testing dataset, which can be used as a green and nondestructive method to timely reflect wheat aging states.


A novel detection method for wheat aging based on the delayed luminescence
Gong Yue-hong 1,2 , Liu Yu-kun 1* , Gong Zhi-le 3 , Zhong Xiao-yan 1 , Zhao Wei-ting 1 , Li Bing 1 , Ge Hong-yi 4 & Lyu Qiong-shuai 1 Wheat aging plays an important role in assessing storage wheat quality and its subsequent processing purposes.The conventional detection methods for wheat aging are mainly involved in chemical techniques, which are time-consuming as well as waste part of wheat samples for each detection.Although some physical detection methods have obtained gratifying results, it is extremely hard to expand their application fields but to stay in the theory stage.For this reason, a novel nondestructive detection model for wheat aging based on the delayed luminescence (DL) has been proposed in this paper.Specifically, after collecting enough sample data, we first took advantage of certain hyperbolic function to fit DL signal, and then used four parameters of the hyperbolic function to feature the decay trend of the DL signal.Secondly, in order to better feature the DL signal, we extracted other six features together with above four features to form the input feature vector.Finally, as the bidirectional long short-term memory (Bi-LSTM) network lacked error-correcting performance, the Bi-LSTM network based on Walsh coding (Walsh-Bi-LSTM) mechanism was proposed to establish the detection model, which made the detection model have the error-correcting performance by reasonably splitting the multi-classification target task.Shown by experimental results, the newly proposed wheat aging detection model is able to achieve 94.00% accuracy in the testing dataset, which can be used as a green and nondestructive method to timely reflect wheat aging states.
Wheat plays an important role in strategic grain reserve in the world by virtue of its excellent storage characteristics.Wheat kernels, as a type of organism, will still carry on respiration function to maintain their lives after harvest 1 .With time elapsing, wheat kernels continue consuming their own energy and then cause wheat aging, which not only deteriorates the quality, but also wastes the quantity of storage wheat 2 .
The main physiological changes of wheat aging were characterized by losing various enzymatic activities and destroying membrane integrality, impairing RNA and protein synthesis, and so forth 3 .Besides the physiological changes mentioned above, it was illustrated that the germination rate of wheat kernels drastically decreased, while the fatty acid value and malonaldehyde content obviously increased after artificial aging 4 .It was shown that the viscosity value and the falling value of wheat kernels would be increased with prolongation of the storage time, which were two disadvantage factors for wheat storage 5,6 .Therefore, the study on wheat aging degree detection is of great importance for ensuring the quality and the quantity of storage wheat.
Generally, conventional detecting techniques for wheat aging mainly rely on chemical methods, such as guaiacol method, triphenyl tetrazolium chloride method, and acidity indicator method.Specifically, guaiacol method is an indirect detection technique to reflect the wheat aging degree via measuring the characteristics of peroxidase which widely exists in animals, plants, and microbes.Wang et al. 7 realized the paddy freshness degree detection using guaiacol method.The principle of triphenyl tetrazolium chloride method is that the embryo colour of vigorous wheat kernels can be dyed to red after certain reactive process provided that dehydrogenase of inner embryo cell absorbs the water thoroughly.Yang et al. 8 successfully applied triphenyl tetrazolium chloride method into detecting wheat aging degree.Acidity indicator method refers to detecting the wheat aging degree based on the regularity of fatty acid value of wheat kernels.Based on the acidity indicator method, Zhao et al. 9 used Methyl Red and Bromothymol Blue as an acid indicator to show the relationship between paddy aging OPEN 1 School of Software, Pingdingshan University, Pingdingshan 467000, China. 2 Henan International Joint Laboratory for Multidimensional Topology and Carcinogenic Characteristics Analysis of Atmospheric Particulate Matter PM2.5, Pingdingshan 467000, China. 3School of Computer and Software, Pingdingshan Polytechnic College, Pingdingshan 467000, China. 4School of Information Science and Engineering, Henan University of Technology, Zhengzhou 450001, China.* email: txlyk@126.com

Sample preparation and instrument parameter configuration
Above all, all the following experiments on wheat samples were carried out according to the institutional guidelines and legislation.Wheat samples in five different years from 2015 to 2019 (we successively use SW2015-SW2019 to label them in the following chapters) were offered by Suiping Grain Bureau, Zhumadian City, Henan province, China.Firstly, we picked out cracking kernels and other foreign materials, and then washed them three times by distilled water.Secondly, we took advantage of the electric blast drying oven to dry the samples till the moisture was (12.5 ± 0.2)%.Thirdly, each type of wheat was prepared about three kilograms and put these five types wheat into five valve bags marked by corresponding year for the subsequent measurements.
The type of UWL analyzer is BPCL-2-ZL, manufactured by Beijing Jianxin Lituo Technology CO., LTD, which can be seen in Fig. 1.The UWL analyzer consists of three parts: ① dark chamber, where we put the tested wheat samples and fix a cold white Light-Emitting Diode (LED) as the induced illuminant; ② photomultiplier (PMT) housing, which is used for changing photon signals into electrical signals; ③ high voltage (HV) supply and pulse counter, taking advantage of the pulse amplitude screening technology and the digital counting function to collect the photon signals.
The type of Digital Lux Meter is AS803, manufactured by Hong Kong Smart Instrument Group CO., LTD.Testing range: 1 LUX ~ 200,000 LUX; operating temperature: − 10-60 °C; measurement error: ± 2%; resolution: 1 LUX.LUX stands for the unite of intensity of illumination.Figure 1 shows the main instruments used in this work.
Compared with other induced illuminants, LED is finally chosen in this work because of the following advantages: • LED belongs to a type of all-solid, anti-seismic, long service life, small single power cold illuminant, and can still work even under low voltage conditions.• LED shows better monochromatic and steady optical properties, and the light intensity is not easily affected by voltage fluctuations.• LED has fast response characteristics, and can timely achieve continuous work states without preheating.
• The volume of LED is small, which is easy for further design and installing in a small place like the dark chamber.

DL signal acquisition
First of all, before carrying out the experiment, keeping the same ambient conditions, such as, indoor temperature at (20 ± 1) °C, relative humidity at (21 ± 6)%, testing time (8:00 a.m.-6:30 p.m.) and so on, is of necessity for the sake of minimizing the influences caused by environmental factors.In addition, we randomly select certain wheat samples to test the analyzer so as to choose the optimum working parameters, and setting the following parameters after many tests: measuring temperature is 20 °C; high voltage is 1030 V; sampling interval is 0.1 s; and the total measuring time for each wheat sample is 100 s.Finally, the empty chamber has been measured several times according to the above-mention parameters in order to obtain the average background noise of the analyzer.
Subsequently, we have installed the LED as the induced illuminant to measure the DL signals of wheat samples in five different years, the schematic diagram of which can be seen in Fig. 2. Specifically, taking wheat sample in year of 2015 for an example, the weight of each sample is (25 ± 0.02) g and the total number of samples is 100 groups, and then we put these samples into a dark box for 30 min before testing in order to minimize the influences caused by other stray light.Furthermore, we have designed a set of button switch device to control the lighting time manually.After many attempts, we have selected the lighting time as 60 s, and then measured DL signals of wheat samples in other years under the same working conditions.

DL signal data preprocessing process
While measuring DL signals of wheat samples, the photon signals are not only influenced by the dark current from other electronic components, but also affected by the changes of environmental conditions.For that reason, the preprocessing to the DL signals is of necessity.We take certain SW2016 sample as an example to show the specific preprocessing process to the DL signal data, which can be seen in Fig. 3. www.nature.com/scientificreports/First, the average background noise intensity (9 counts/0.1 s) of the UWL analyzer is calculated according to the experimental parameters mentioned above, which is seen in Fig. 3a.Second, the signals are collected according to the mode of subtracting the average background noise intensity value, as shown in Fig. 3b.Third, the DL signals obtained in the second step are smoothed and denoised by the sliding average method to get the final sample data, and the sliding window value is set to 5, which is shown in Fig. 3c.

Classification feature extraction method
In order to better feature DL signals of SW2015-SW2019, we take advantage of following methods to extract ten classification features in this work.

Hyperbolic relaxation feature extraction method
It has been validated that the organism system induced by a white or other monochromatic illuminants shows a dynamic relaxation behavior, which has the following characteristics: • show a comparative long relaxation time.
• the end of relaxation process is proximate to SL intensity.
• the whole relaxation process can be fitted by a hyperbolic function perfectly.
The detailed hyperbolic function can be described as where I 0 is the initial luminescence intensity, which depends on the properties of sample as well as lighting condi- tions; τ is a time feature, which is only related to the properties of sample; β used for controlling the relaxation rate is an exponential factor 22 .Figure 4 shows a schematic plots of DL signal of certain sample SW2016 fitted by hyperbolic function.
Table 1 shows the hyperbolic fitting parameters and evaluation indices of the DL signals of certain samples in five years, where SSE stands for Sum of Squares due to Error, and RMSE is the abbreviation of Root Mean Squared Error, and R 2 represents Coefficient of Determination.Seen from the evaluation indices in Table 1, the value of RMSE is less than 21 and the R 2 value is more than 0.99, which validates that the DL signals of wheat samples can be perfectly fitted by certain hyperbolic functions, and lays a foundation for feature extraction to DL signals based on the fitting functions.
Moreover, we can also resort to a time window to calculate the integral calculus of (1) with the range of [0,T], which is shown as in formula ( 2) where, T is the measuring time.E(T) stands for the calculus intensity of the DL with the range of [0, T] , and the value of E(T) is depended on the parameters I 0 , τ , β and T. (1)

Other feature extraction methods
Besides the above-mentioned four features, we also elaborately choose the following six features to shape the final classification feature vector to feature the DL signals.

Instantaneous frequency
Instantaneous frequency (IF) of a time sequence signal is that the mean of frequencies f in the signal changes with diverse time instant parameters t , which provides a measure of signal energy concentration in the frequency domain as a time function 23,24 .Thus the function IF(t) can be used to estimate the IF of a signal at instant time t through computing (3)   where P(t, f ) is the signal power spectrum.
In this work, due to the DL signal is measured in a time interval 0 ≤ t ≤ T rather than in an infinite length.For this reason, we take advantage of a special method for computing the power spectrum of time sequence x k , k = 0, 1, ..., N − 1 , and then the power spectral density of x k can be approximately described as 25 .where t stands for the sampling interval.
We can rewrite (4), assuming f = j f , �f = 1/(N�t) , t = "1", and then we can obtain where X k is the discrete Fourier transform (DFT) of x k , which can be seen as follows. (3) x k e −2πifk�t

2
(5)  In order to fully use the power spectrum characteristics, we resort to the Kaiser window to scale P j .Thus (5) can be rewritten as Among which where w j is the weight of the Kaiser window function.
So the P j in ( 7) can be calculated by the following steps: • In order to simplify the computation process, the length of initial time series is taken N = 2 n , if not, padding that with zeros, n stands for a positive integer.• Use a corresponding window function w k to weight the above series.
• Resort to the fast Fourier transform (FFT) to calculate the discrete Fourier transformation (DFT) of the weighted series w k x k .• Calculate P j by (7).

Spectral entropy
Spectral entropy (SE) is used to measure the signal spectral power distribution 26 .Normally, there remain two steps for calculating the SE at time t , which is denoted as H(t) .Firstly, assume the probability distribution of a time-frequency power spectrum as P(t, f ) , and p(t, m) stands for frequency point m, m = 1, 2, ..., N , which can be calculated by where f ∈ [0, f s /2] , f s denotes the sampling frequency, and we take f s = 10 HZ in this work due to the sampling time is 0.1 s.Secondly, according to the computation process of Shannon entropy, the spectral entropy H(t) can be obtained by

Approximate entropy
The approximate entropy (ApEn) algorithm was initially proposed by Pincus, aiming at measuring the properties of random series 27 .Two striking advantages of the ApEn lies in its lower dependency on the length of initial time series and robust resistance to the noise carried by the initial data.The detailed computing process of the ApEn algorithm is: 1. Partition the initial series X = {x(i), i = 1, 2, ..., N} into an m-dimensional vector u(i) , which is described as Here, m stands for the dimension of the pattern vector, and N is the initial length of the time series.2. Calculate the distance d[u(i), u(j)] between vector u(i) and vector u(j) using (12).
3. Count the numbers of d[u(i), u(j)] < r , where r is a positive real number, which is known as the similar tol- erance threshold value.Consequently, calculate the proportion of d[u(i), u(j)] < r among the total number of vectors, which is labeled as C m i (r) 4. Calculate the logarithm of C m i (r) , and then, obtain its mean value by (14).Here, the mean value is labeled as H m (r).x k e −2πi jk N , j = 0, 1, ..., N − 1 www.nature.com/scientificreports/ 5.By increasing the dimension from m to m + 1 and repeating above steps 2-4, we can obtain H m+1 (r).
Then the definition of the ApEn can be described as If N is finite, formula (15) can be rewritten as follows: From the above formulas, we can deduce that the ApEn value is prone to be lager provided that the initial series is more complex.Thus, the ApEn can be used as a classification feature to characterize the DL signals of wheat samples.

Median
Rearranging the measured DL data of a sample in a descend order, the value in the middle of the arranged sequence is the median.Namely, in the entire sequence, the number of less than the median and greater than the median is equivalent.If the total number is even, take the average of the most middle two values as the median.

Quartile deviation
Quartile is a set of data in a descend number, using three points to divide the whole data into four equal parts.We mark the corresponding three division points by Q1,Q2 and Q3 respectively.Q1 (the first quartile) indicates 25% of the data less than or equal to Q1; Q2 (the second quartile, namely median) means 50% of the data is less than or equal to Q2; Q3 (the third quartile) points 75% of the data is less than or equal to Q3.Among them, the distance between Q3 and Q1 is also called quartile deviation (QD), denoted by QD = Q3-Q1.The QD reflects the dispersion degree of the middle 50% data, and the smaller of QD reflects the more concentrated the data in the middle is, and vice versa.Sometimes, we will resort to interpolation method to calculate the value of Q1 and Q3, provided that the number of (n + 1) can not be divisible by four, where n is the total sampling points.

Mean deviation Mean deviation (MD) can be calculated by the following equation:
Here N is the total sampling points, and µ is the corresponding mathematical expectation.The value of MD characterizes the differences degree between each sample point and the arithmetic mean.The greater value shows the smaller representation of the arithmetic mean and vice versa.

Bi-LSTM network
Long Short-term Memory (LSTM), proposed by Hochreiter and Schmidhuber 28 , is a special Recurrent Neural Networks (RNN), which is capable of handling long-term dependencies as well as solving the gradient vanishing problem existed in conventional RNNs.To the LSTM network, it takes advantage of diverse gates shaped by input data multiply different vectors to control data flows.Concretely, the input data controls how much the new data can be used in the current memory cell; meanwhile, the forget gate makes a decision that how much data should be removed from the former memory cell; finally, the output gate decides how much data will be output from the current memory cell.The detailed inner structure of LSTM units is shown in Fig. 5, and the red activation function stands for hyperbolic tangent function while the yellow activation function is sigmoid function.By virtue of these diverse function gates, LSTM network is capable of processing time sequences effectively to achieve the purpose of classification or prediction.
To initial DL data, we resort to a zero mean normalization method (Z-score standardization) to normalize the initial data.The normalization equation can be described as where x i (1 ≤ i ≤ t) stands for all the elements of the input vector, µ= 1 t t i=1 x i and σ = represent its corresponding mean and sample standard deviation respectively.
Consequently, we take advantage of LSTM network to train the labeled normalization DL data x = (x 1 , x2 , ..., xt ) , and the whole process can be described by the following formulas according to Fig. 5 (14) where W * stands for the corresponding matrix weights;b * denotes the different bias scalars; σ (•) is sigmoid function; tanh(•) represents hyperbolic tangent function.i, f , o, g, s are input gate, forget gate, output gate, can- didate values, and new cell state successively; y t refers to the final output and ⊙ means the element-wise product between two vectors.
Softmax Classifier: the softmax function s = [s 1 ,s 2 ,...,s M ] can be defined as (25), which maps M input data vector into M normalized output data.
where w i stands for the weight vector of the fully connected layer, k f refers to the final output in the second layer, (•) T means transpose operation.
In order to train the weights of LSTM network and measure the difference between the normalized output and the true label, we have introduced cross entropy as the loss function L(w) , and adopted L2 regularization method so as to prevent over-fitting and improve the generalization performance of network, which can be describes as where y i denotes for the true labeled data for the concrete wheat storage year, and w i is the weight parameter, and ∈ R represents the regularization coefficient in order to adjust weight between regularization term and original loss value.Meanwhile, M stands for the categories of output of softmax classifier, and N is the number of training set.
To the error term of LSTM, it has been involved in two directions: one way is that we can compute it along the time dimension by back propagation, namely, computing every error term before t moment; the other way is transmitting the error to the front layer.However, fewer researchers pay attention the second one, so we mainly resort to the first method to calculate the error term, which can be described as where δ l i,t , δ l f ,t , δ l g,t and δ l o,t are standing for the error term of input gate, forget gate, candidate value and output gate in l th layer respectively.Furthermore, we can deduce them by following formulas.Here, f ′ (•) means the differential of active function, and we always choose tanh and sigmoid function as active function in LSTM network.Thus, each element in weight matrix is less than one, and we can see from the above formula that there remain positive or negative signs among the four terms, which makes the value of δ l−1 t decrease gradually and causes the problem of gradient vanishment.
After taking all the above-mentioned factors into consideration, we finally resorted to the Back Propagation Through Time (BPTT) algorithm to train the parameters of LSTM network based on the above loss function 29 , in which Adam Optimizer was used for the sake of making the computation process more efficient 30 .
In order to improve the classification performance of the proposed model, we have introduced a bidirectional LSTM (Bi-LSTM) network to learn the "future information" of input sequence 31 .Figure 6 shows the concrete Bi-LSTM network structure, among which {x 0 , x 1 , x 2 , ..., x t } and {y 0 , y 1 , y 2 , ..., y t } are the input and the output sequence respectively; A and Á stand for the respective LSTM units; s 0 ,s t ,s 0 ′ and s t ′ are the different cell states.This type of network not only can be trained in two directions simultaneously, but also has the separate hidden layers.

Mechanism of encoder-decoder based on Walsh code
In order to make the classification process more robust, we resort to Walsh code to generate the labeled output sequences in this paper.There exist several methods for producing Walsh code sequence 32 , and we mainly take advantage of Hadamard matrix to generate Walsh code in this work.The concrete steps can be seen as follows: we assume H N denotes the Hadamard matrix with dimension N × N , among which any two rows or columns are orthogonal.Namely, if we regard the row or column in H N as a function, the values of cross correlative func- tion between them are zeros.Walsh code can be generated by the following recursive procedure, which can be describes as Here = 2 n , n ∈ N + ,H N is the logical reverse of H N , and normally we set H 1 = [0][0].Since our final classification targets in this paper are five, we take N = 8 to generate the target sequences, which can be described as We can obviously see from (33), either the ith row or the ith column ( i is row or column number of the matrix) in H 8 shows the same coding pattern, which not only exhibits a strong robust ability, but also easily extends according to the requirement.In this work, we successively select the first five rows or columns as the representations of the corresponding wheat sample labels from 2015 to 2019.

Walsh-Bi-LSTM network
As the Bi-LSTM network lacks error-correcting performance, we propose the Walsh-Bi-LSTM network to establish the final detection model, aiming at making the detection model have the error-correcting performance by reasonably splitting the multi-classification target task.The whole flowchart of the detection model is shown in Fig. 7.Moreover, in order to test the error-correction performance of the model, noise is added to all data before entering the model, where the dashed box represented the Walsh encoding module.

Experiment results
In this part, we carried out comparable experiments to validate the classification performance of our proposed model.In order to make a comparison, we tested the wheat aging classification performance by the normal Bi-LSTM network and the Walsh-Bi-LSTM respectively.The performance of classification model was assessed by three indices: training accuracy curve, training loss curve and confusion chart.

Datasets
We carried out all the experiments based on the DL data of wheat samples stored in five different years.First of all, we take advantage of above-mentioned feature extraction methods to obtain the final datasets.In addition, the datasets of each type were divided into two groups, namely, training group (80 samples) and testing group (20 samples), and each sample was marked by the corresponding classification labels.

Implementation details
In this work, we adopted normal Bi-LSTM network and Walsh-Bi-LSTM as the detection models respectively for the sake of making a comparison between their classification performance.After several tests, we finally chose the following hyperparameters: The batch size was 8; the number of epochs was 100; learning rate was 5e−4; hidden layers was 128 and Adam optimizer was adopted.To the Walsh-Bi-LSTM, we took H 8 to encode and decode the output sequences, among which we chose the top five rows or columns as the representation of corresponding five classification labels.Furthermore, we used the ten above-mentioned features extracted from DL signal as the training and testing data for two networks.Subsequently, both of networks was implemented by TensorFlow framework and trained on the basis of NVIDIA GeForce GTX 1060 3 GB.

Classification evaluation of Bi-LSTM network and Walsh-Bi-LSTM network
In this part, we have trained the Bi-LSTM network and the Walsh-Bi-LSTM network separately and made a classification performance evaluation between them.Seen from Fig. 7, we have taken z-score operations and added noise so as to test the robust performance of classification models.After simulating by computer, the training average accuracy of the Bi-LSTM network and the Walsh-Bi-LSTM network is 87.3% and 99.6% respectively, which can be seen in Figs. 8 and 9. To the testing dataset, the average accurate classification rate of the Bi-LSTM is 87.0%, and the misclassification results mainly scatter in the adjacent years due to the one-hot vector used as final output labels, which can be seen in Fig. 10.Meanwhile, the average accurate classification rate of the Walsh-Bi-LSTM in the same testing dataset is 94.0%, and the misclassification results randomly distribute in other years due to the Walsh encoding mechanism, which can be seen in Fig. 11.
Experimental results show that the average training classification accuracy and testing classification accuracy of Walsh-Bi-LSTM network achieve 99.6% and 94% respectively, which have been improved by 12.3% and 7.0%    www.nature.com/scientificreports/Seen from Table 2, we observe that the Bi-LSTM network can obtain a better classification results both in training set and testing set, and the main reason is that DL signals of wheat kernels have strong time association features.Moreover, it is also validated that Walsh-Bi-LSTM network proposed in this work does have errorcorrecting performance, which is a feasible and reliable way to detect wheat aging degree.

Conclusions
We have proposed a set multi-class wheat aging detection model using the Walsh-Bi-LSTM network based on DL signals of wheat kernels.Shown by the experiment results, the DL signals of wheat kernels can be used as an index to reflect their inner physiological activity like wheat aging.The detection technique to wheat aging based on DL is an eco-friendly and nondestructive method.Moreover, compared with detection methods based on SL signals, the detection methods based on DL signals have a higher SNR as well as easily extract classification features.Since DL signal obeys certain hyperbolic function decay trend, we can resort to Bi-LSTM network and Walsh-Bi-LSTM network to establish the classification model respectively after obtaining enough data samples.The simulation results show that the Walsh-Bi-LSTM network is prone to be more efficient and feasible in detecting wheat aging degree.Our next work is establishing a model for predicting the wheat aging states, which can help granary managers to evaluate the stored wheat aging states and make scientific decisions.

Figure 1 .
Figure 1.The main instruments used in this work.

Figure 2 .
Figure 2. The schematic diagram of measuring DL signal of wheat samples.

Figure 3 .
Figure 3. Data preprocessing process to the DL signal of SW2016 wheat sample.

Figure 4 .
Figure 4. Schematic plots of DL signal of certain sample SW2016 fitted by hyperbolic function.

Figure 5 .
Figure 5.The detailed inner structure of LSTM unites.

Figure 7 .
Figure 7.The whole flowchart of the detection model for wheat aging.

Table 1 .
The corresponding fitting parameters and evaluating indices of DL signals in five years.

Table 2 .
The classification results of different models based on the same features.